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ABSTRACT 

Binaries containing a stellar-mass black hole and a recycled radio pulsar have so far eluded detection. We 
present a focused investigation of the formation and evolution of these systems in the Galactic disk, highlighting 
the factors that limit their numbers and the reasons why they may be extremely rare. We surmise that the 
birthrate of black-hole/recycled-pulsar binaries in the Galactic disk is probably no higher than ~10 _7 yr _I , and 
may be much less, including zero. Simple arguments regarding common-envelope evolution suggest that these 
binaries should have orbital periods less than lOhr and an average lifetime of <10 8 yr before coalescence due 
to the emission of gravitational radiation. We expect that fewer than ~10 of these compact, relativistic binaries 
currently reside in the Galactic disk, less than 0.1-1% of the number of double neutron stars. The discovery of 
two or more black-hole/recycled-pulsar binaries using current radio telescopes would tightly constrain certain 
ideas regarding the evolution of massive stars, dynamical mass transfer, and black-hole formation. 
Subject headings: binaries: close — black hole physics — gravitational waves — pulsars: general — stars: 
neutron 



1. INTRODUCTION 

Thirty y ears have passed sinc e the discovery of PSR 
B 19 13+ 16 JHulse & Tavlorlll975l) . which has the dual dis- 
tinction of being the first known binary radio pulsar and dou- 
ble neutron star. The intervening years — particularly the last 
decade — have seen enormous growth in the observed number 
and variety of binary radio pulsars. Among the ^100 known 
systems (see Stairs 2004 for a representative listing and ref- 
erences), orbital periods range from 90 min to several years, 
most have white-dwarf companions with masses from sub- 
stellar to near the Chandrasekhar limit, three have a stellar 
companion, and eight probably contain a second neutron star 
( Stairs l2004 lFaulkner et al.l2005l) . 

With a handful of exceptions, most radio pulsars in binary 
systems are recycled, characterized by low inferred surface 
magnetic field strengths of <10 10 G and spin periods of 2- 
100 ms. The term "recycled" refers to the accretion of matter 
and angular momentum by the neu tron star, whereby the spin 
period is reduced to small values (|Smarr & B landfordl ll976t 
lAlpar et al]l!982t ljoss"& RaDParjortl ll983l) . In contrast, the 
majority of isolated pulsars (hereafter called "normal") have 
magnetic fields of 10 n -10 I3 G and spin periods of 0.1-1 s. 
Recycled pulsars show much greater stability as clocks, al- 
lowing for very precise dynamical measurements if the pulsar 
is in a binary. Precision timing of a recycled binary pulsar 
quickly reveals its orbital motion, and if the orbit is suffi- 
ciently compact or inclined, relativistic gravity may cause sig- 
nificant deviations from a Keplerian model, as in the classic 
example of PSR B1913+16. 

The eight observed double neutron stars (DNSs) have a spe- 
cial status among binary pulsars. As exotic endpoints of mas- 
sive binary stellar evolution, they probe important details of 
stellar physics. Moreover, since the bright pulsar is recycled, 
and both components are essentially point masses, a subset of 
these systems have proved to be extremely rich laboratories 
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for studying post-Newtonian gravity. The newly discovered 
double pulsarbinary, PSR J0737-3039, stands out as having a 
nearly edge-on orbit a nd an orbital period of 2.4 hr, the short - 
est among the DNSs llBurgav et al J 120031 iLvne et alJl2004l) . 
Finally, the known DNSs, especially PSR J0737-3039, im- 
ply a large number of similar systems in the Milky Way and 
distant galaxies that will produce in the next ~1 Gyr a power- 
ful burst of electromagnetic and gravitational radiation in the 
final throes of coalescence (e.g., Kalogera et al. 2004). 

All the factors that make DNSs exceptional are greatly en- 
hanced if the pulsar's companion is replaced with a canonical 
10M Q black hole. To date, there has been no detection of 
a binary containing a stellar-mass black hole and a recycled 
pulsar (hereafter BHRP). Several earlier theoretical studies 
provi de insight into the formation of B HRPs in the Galactic 
disk. Naravan, Piran, & Shemi (1991) first mentioned these 
systems, identified a likely formation scenario (see § [2j, and 
suggested a birthrate of <10~ 6 yr based on t heir non detec- 
tion. B inary population syn t hesis s t udies bv iLipunov et alJ 
i 1994, ISipior & Sigurds"sor] (120021), IVoss & Taurisl (12003 ), 
and Sipior. Portegies Zwart. & Nelemans (2004) theoretically 
constrain the BHRP birthrate and the expected orbital pa- 
rameters. Each of these papers also addresses the forma- 
tion of systems where the radio pulsar is not recycled (see 
also [Belczvnski et al.l l2002l) . which are not considered here. 
These latter systems are certainly interesting, and probably 
much more numerous than BHRPs, although they they are not 
nearly as desirable as astrophysical laboratories. We also do 
not investigate the formation of BHRPs in globular clusters 
(Sigu rdssonl2003l) . 

Our main goals are to (i) dissect the most important steps 
and uncertainties in BHRP formation (§§ [3 and [3J, (ii) con- 
strain the BHRP birthrate in the Galactic disk (§0), (in) in- 
dicate the range of binary parameters for newly formed sys- 
tems and the current-epoch population (§ |5j, and (iv) eval- 
uate the likelihood of discovering BHRPs in pulsar surveys 
(§ EJl. Each of these topics has received some attention in 
one or more of the above papers, but no clear picture has yet 
emerged, and previous work has not examined all the major 
physical uncertainties. We utilize simple, semi-analytic calcu- 
lations to present a complete portrait of BHRP formation. The 
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next section gives a preview of the major evolutionary hurdles 
in BHRP formation, and is followed by a more detailed and 
quantitative discussion in §|3] 

2. OVERVIEW OF FORMATION AND UNCERTAINTIES 

In this paper, we consider the most obvious formation 
channel for BHRPs in the Galactic disk, which is qual- 
itatively identical to the standard model for forming the 
known DNSs (e.g.. iBh attacharva & van den Heuvel| 1199 It 
Tauris & van den Heuvelll2003l) 4 . The scenario begins with 
a massive primordial binary, where the primary star is the 
neutron-star progenitor of mass ~8M Q to >25M Q . Dur- 
ing a phase of dynamically stable mass transfer from the 
primary, the initially less massive secondary reaches a mass 
of >25Mq, and will ultimately leave a black-hole remnant. 
Mass transfer ends when the core of the primary is exposed. 
The core then explodes as a Type Ib/c supernova, leaving be- 
hind a neutron star that may have received a large impulsive 
kick. Before the secondary evolves significantly, the neutron 
star may be a radio pulsar or an X-ray source accreting from 
the wind of its massive companion. Because of the extreme 
binary mass ratio, the second episode of Roche-lobe over- 
flow is dynamically unstable. The neutron star is engulfed 
by the stellar envelope and dragged toward the secondary's 
core. Following the successful expulsion of the common en- 
velope, the neutron star orbits the core, which then collapses 
to a black hole; for less massive secondaries, a second neutron 
star is formed at this stage. Either before, during, or after the 
spiral-in, the neutron star accretes a small amount of mass and 
is recycled. We now address several fundamental uncertain- 
ties in the evolution of the massive secondary, spiral-in phase, 
recycling process, and formation of the black hole. 

2.1. Stellar Evolution Uncertainties 

Surprisingly, the evolution of massive stars is still quite 
poorly understood. It is generally estimated that the min- 
imum mass of a single star that ends its life as a black 
hole is ~25M Q (see Podsiadlowski, Rappaport, & Han 2003 
for discussion and references; hereafter PRH03). However, 
iBrown et al.l (|2001) show that when a star less massive than 
40-50M Q loses its hydrogen-rich envelope before helium 
core burning (case B mass transfer), its core evolution is 
changed dramatically, and the star will most likely leave a 
neutron-star remnant. Therefore, mass transfer from the sec- 
ondary in the above BHRP formation scenario must begin 
after it has completed most or all of it s helium core burn- 
ing ph ase (case C mass transfer; see also lWellstein & Langerl 
1999). The formation of BHRPs is then very sensitive 
to the radial evolution of the massive s econdary (see also 
iPortegies Zwart. Verbunt. & Ergmal [T997h and uncertainties 
in the stellar physics. 

The range of binary separations for case C mass transfer 
depends mainly on the difference between the maximum ra- 
dius of a single star and the radius at the tip of the first giant 
branch. For massive stars, this range may be as small as a 
few percent of the radius in some evolutionary models or a 
factor of >2 in others. The result depends on when stars un- 
dergo helium core burning (see, e.g., the difference between 
the evolutionary tracks in Fig. 1 and 2 of Pfahl, Rappaport, 

4 The only alternative scenario for DNS formation, involving the spiral- 
in of two stellar cores inside a common envelope (Brown (1995), is proba- 
bly not a viable channel for BHRP formation, since the pulsar would be the 
second-born compact object and hence not be recycled (see also J. Dewi, Ph. 
Podsiadlowski, & A. Sena 2005, in preparation). 



& Podsiadlowski 2002; hereafter PRP02). If helium burning 
occurs only while the star is a blue supergiant, the range of 
case C mass transfer is greatly increased compared to when 
the star burns helium as a red supergiant (by perhaps a factor 
of 4 or more; Dewi et al. 2005, in preparation). The ob- 
served distribut ion of supergiant s in the Hertzsprung-Russell 
diagram (e.g.. iHumphrev s 1984) suggests that stars as mas- 
sive as ^40 Mq spend a significant fraction of their helium 
burning phase as both blue and red supergiants. None of the 
present theoreti cal models is yet fully co nsistent with the ob- 
servations (e.g.. lLan ger & Maeder 1995). 

2.2. Severity and Implications of the Spiral-In 

For a neutron star sinking into the envelope of a >25M Q 
star, the orbit may have to shrink by a factor of >10 3 before 
the envelope can be expelled, according to recent calculations 
of envelope binding energies (Dewi & Tauris 2001; PRH03). 
In order to avoid a merger of the neutron star and stellar 
core of radius ~17? Q , the initial orbital separation must be 
>2OOO/?0. Stars of m ass 25-50Mm attain maxim um radii of 
~1500-2500/? Q (e.g.. lHurlev. Pols & Toul2000l) . and so the 
range of orbital separations for successful systems is restricted 
to ~2000-4000fl Q . A star that fills its Roche lobe in such an 
orbit is very evolved, and in most cases has passed through 
the ph ase of helium core burning. Thus, in the IBrown et al.l 
(2001) picture, we are safe in using ~25M Q as the mass 
threshold for black-hole formation. 

If the orbits contract by a factor of 1000, the final separa- 
tions are <ARq. However, if the contraction factors are gen- 
erally >2000, which is quite plausible, then it may be that 
no BHRPs are formed. Another serious complication arises 
because the secondaries of interest are very evolved prior to 
the spiral-in. Prodigious stellar winds during the late phases 
of evolution may generally cause the secondary's Roche lobe 
to expand fas ter t han the star, precluding mass transfer alto- 
gether (see § 13.41 Portegies Zwart, Verbunt, & Ergma 1997; 
Kalogera & Webbink 1998; PRH03). On the other hand, if the 
star can catch its Roche lobe after losing a significant amount 
of mass, the orbit is wider at the start of the spiral-in and the 
envelope mass that must be ejected is reduced, giving the sys- 
tem a better chance to survive.. 

We assume that the neutron star does not accrete at hyper- 
critica l rates and collapse t o a black hole during the spiral-in 
(e.g.. lChevalierll993lll99d) . If this were to occur generally, as 
may, in fact, be the case, then neither DNSs nor BHRPs would 
form by the channel being discussed (e.g., Brown 1995). 

2.3. Recycling 

The neutron star in a BHRP is mildly recycled after accret- 
ing a small amount of matter (<0.001M Q ) from, e.g., the 
wind of its companion prior to the spiral-in, the stellar en- 
velope during the spiral-in itself, or a residual accretion disk 
after the orbital contraction has halted. For most DNSs, recy- 
cling may result during a short phase of stable mass transfer 
from the Roche-lobe filling core of the secondary if it is suffi- 
ciently l ow in mass and if it is exposed before ce ntral helium 
ignition JPewi & Polsll2003l llvanova et alJ l200l . However, 
this process does not apply to BHRPs, because secondary's 
core is too massive for the transfer to be stable. 

An interesting possibility is if recycling in BHRPs occurs 
mainly after the spiral-in, as the neutron star accretes a frac- 
tion of the strong wind from the secondary's core. Since the 
core has most likely passed through the phase of central he- 
lium burning, only <10 4 yr remains before it collapses. At the 
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Eddington-limited accretion rate of ~4 x 10~ 8 Moyr _1 (for 
pure helium), the neutron star can accept only ^10~ 4 Mq of 
material. The short timescale, low accreted mass, and possi- 
bly low angular momentum of material accreted from a wind, 
all suggest that pulsars in BHRPs may not be strongly recy- 
cled, and could have systematically longer spin periods and 
higher surface magnetic fields than pulsars in DNSs (see Ar- 
zoumanian, Cordes, & Wasserman 1999 for a related discus- 
sion and references). 

2.4. Collapse of the Secondary 

Uncertainties in black-hole formation include the effects of 
winds from the massive progenitor and its exposed core in a 
binary system, stellar rotation, and the dynamics of core col- 
lapse (e.g., Fryer & Kalogera 2001 and references therein). 
For the collapse, there are two possibilities. Either the core 
collapses promptly with out being accompanie d by a bright 
supernova (e.g.. lGourgoulhon & Haensell 19931) . or a neutron 
star forms first and subseque ntly collapses to a black ho le after 
accreting fallback material (Wooslev & Weaver 1995), eject- 
ing a significant fraction of the envelope and resulting in a 
successful supernova. The two-stage scenario has been in- 
ferred for the black-hole progenito r in the binary Nova Scor- 
pii (GRO J 1 655-40: iPodsiadlowsk i et al.ll2002l) . Moreover, 
the high space velocity of Nova Sco requires that the black 
hole received a large natal kick of >100kms _1 (Brandt, Pod- 
siadlowski, & Sigurdsson 1995). It is not k nown if most black 
holes receive significant kicks a t birth (e.g.. lFrver & Ka logera 
l2001tUonker & Nelemansl2004l) . A plausible guess is that the 
kick momentum received by newbor n bla ck holes is similar to 
that received by neutron stars (see § !3.6t . 

3. FORMATION DETAILS 

Here we elaborate on the points in the last section. Each 
stage in the formation of a BHRP is described explicitly 
and quantitatively, although with a more transparent approach 
compared to many previous binary population synthesis stud- 
ies. Our chosen assumptions and prescriptions for stellar evo- 
lution are fairly standard, but are not unique. Further details 
and references may be found in PRP02. 

3.1. Primordial Binaries 

In our standard model, the mass of the primary is drawn 
from a power-law initial mass function, p{M™) oc (M™)~ 2 5 , 

where M™ is in the range of 8M Q to M t h, the mass threshold 
for black-hole formation. Because of the uncertainty i n M t h 
and its possible dependence on binary evolution (see § 12. 1> . 
we consider two scenarios, one where M t h = 25 M Q (scenario 
I), and the other where M± = 40M© (scenario II). The sec- 
ondary mass M™ is taken from a flat distribution in mass ra- 
tios, q PB =M™/M™ < 1. Alternative choices for/?(M™) and 
piqps) are considered in §|5]and Table 1. Without much loss 
in generality, we assume that the orbit is circular, with a sepa- 
ration, ape, drawn from a logarithmically flat distribution that 
extends from ^10/?© to a somewhat arbitrary maximum of 
10 6 R Q . 

3.2. Stable Mass Transfer 

In isolation, the primary would grow to a maximum radius 
of > 1000,^0 before exploding as a supernova. Mass transfer 
from the primary o ccurs if it fills its Roche lobe of radius = 
apBrhiqpB), where (Eggleton 1983) 

0-49 (1) 
0.6 + ^111(1 +^b /3 ) ' 




FIG. 1. — Cumulative distribution of secondary masses following stable 
mass transfer, where the primary initial mass function is p(M™) oc (M™)" 2 5 
for SMq < < M t h, and the secondary mass is drawn from a flat distri- 
bution in mass ratios such that Mf B /Mp B > 0.5. The mass capture fraction is 
fixed at (i = 0.5 (dashed curves) or 1.0 (solid curves). Each pair of curves for 
a given j3 corresponds to a different mass threshold for black-hole formation 
from the primary: = 25 Mq (scenario I) or = 40 Mq (scenario II). 



with typical values of r^ = 0.5 ±0.1. Mass transfer occurs 
over ~2 decades in «pb, comprising c±40% of the binaries. 

Mass transfer may be stable or dynamically unstable de- 
pending on and the evolutionary state of the primary. If 
mass transfer is to be stable, then qps cannot be too small, and 
the primary should not be a deeply convective red supergiant 
(see PRP02 for a discussion, caveats, and references). In this 
regard, we assume that ^pb > 0.5, which is true for half of the 
binaries. By excluding very evolved primaries at the onset of 
mass transfer, we limit the stellar radius to less than ~200- 
1000/? Q , where larger values correspond to higher masses. 
We also neglect systems where the primary transfers matter on 
the main sequence («pb < 30— 40/?q), which will not evolve 
into the objects of interest here. With these cuts, the net frac- 
tion of massive binaries that undergo stable mass transfer after 
the primary has left the main sequence is /mt ~ 0.1. 

The primary is expected to lose its entire envelope during 
mass trans fer, leaving its hydr ogen-exhausted core of mass 
(e.g. JHurlev. Pols. & Toul2000l) 

M^~Q.\{Mf/M Q ) l35 M C3 . (2) 

Let (3 be the fraction of transferred matter accreted by the sec- 
ondary. The final secondary mass is then 

Mf T = M™ + /3(AT™ -M CiP ) . (3) 

We assume that the secondary is reset on the main sequence 
appropriate for its new mass, which is not justified for qp% 
near 1. 

Given (3, it is straightforward to determine the distribution 
of M^ T given the distributions of qpg and M™. Figure ^ 
shows the resulting cumulative distribution of Af, for (3 = 0.5 
and 1 .0 in scenarios I and II. The fraction of systems in which 
the secondary collapses to a black hole and that remain po- 
tential BHRP progenitors can be es timated from Fig. ^ If> 
in accord with the remarks in § 12.21 we assume that all post- 
transfer secondaries of mass >25M leave black-hole rem- 
nants, then for (3 = 0.5-1, the fraction of systems in sce- 
nario I is /bh — 0.04-0.16, while in scenario II the fraction 
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FIG. 2. — Supernova survival probability as a function of a k /v b for a FIG. 3. — Cumulative distribution of circularized, post-supernova orbital 

Maxwellian kick-speed distribution, where a k is the one-dimensional ve- separations, calculated under the assumptions that the pre-supernova binary 

locity dispersion. The probability is <10% for a k /v b > 1.7 and <1% for mass i s 4OM , 20% of the binary mass is lost in the explosion, and the 

a khb > 3.9. For a k /v b > 2, the probability approaches the cubic trend distribution of pre-supernova orbital separations is logarithmically flat from 

~64%(a k /v b r 3 (dashed line). 0.2AUtolOAU. 



is /bh — 0.12-0.23. Setting the minimum secondary mass to 
2OM , we find that /bh spans the range of 0.2-0.37 and 0.13— 
0.3 in scenarios I and II, respectively. Likewise, a minimum 
mass of 30M Q gives fractions of 0.07-0.14 and 0.006-0. 14 in 
the two scenarios. 

The orbital separation following stable mass transfer, gmt, 
depends sensitively on (3 and the specific orbital angular mo- 
mentum taken by escaping matter. A plausible first approxi- 
mation is that flMT ~ flPB- The resulting logflMT distribution 
is then roughly flat from amt ~ 0.1 AU to several AU, with 
a declining tail that extends to ~10AU, which is correlated 
with higher post-transfer secondary masses. For separations 
>10AU, the primary cannot fill its Roche lobe during the 
mass transfer phase. 

3.3. Supernova of the Primary 

After mass transfer, <1 Myr remains before the Type Ib/c 
supernova explosion of the primary's exposed core and for- 
mation of the neutron star, for which we assume a mass of 
Mns = 1 AMq. Radial expansion of the naked core of the pri- 
mary and wind mass loss do not affect our main results and 
are neglected. Associated with the core collapse and super- 
nova is impulsive mass loss and a kick to the newborn neu- 
tron star. Measured proper motions of isolated radio pul- 
sars suggest a median kick speed of lOO-SOOkms" 1 (e.g., 
i Hansen & Phinnevlll997l lArzo umanian . Chernoff. & Cordesl 
2002). The functional form of the underlying kick distribution 
is weakly constrained. Moreover, there are observational and 
theoretical reasons to believe that the average kick magnitude 
may be much smaller for many neutron stars born in binary 
systems (Pf ahl et al] 120021 Ipodsiadlow ski et alJ 12004). We 
now quantify the dynamical effects of supernova mass loss 
and kicks on potential BHRP progenitors with Mf T > 25 M Q . 

Since the primary's core mass is <8M Q , for BHRP progen- 
itors the fractional binary mass lost in the supernova is 



A - Mc ' p 
Ans = 



-M NS 



<0.2 



(4) 



and so the dynamical effect of mass loss is relatively small. 
There is a significant probability for a binary to be unbound 



by the supernova if the neutron-star kick speed, v k , is compa- 
rable to or larger than the relative orbital speed of the binary 
components before the explosion: 



Vb : 



190kms _1 



( M \ 



MT 



1 AU 



V.40M, 



flMT 



1/2 



(5) 



where Mf T = M™ T +M c . p . When v k /v b = 1, the probability 
that the binary remains intact is roughly 50% if the kick di- 
rections are distributed isotropically. To be more quantitative, 
we must adopt a specific distribution. 

A Maxwellian distribution, p(v k ) oc v|exp(-v|/2crf), has 
the attractive features of being a reasonable characterization 
of the pulsar velocity data (for ak = 100-200kms _1 ) and 
having only one free parameter. Figure |2] the result of a 
simple numerical integration, shows the survival probabil- 
ity versus a k /vb for an assumed fractional mass loss Ans = 
0.2. In Fig. |2j the probability asymptotically approaches 
~0.64(cri:/vi)~ 3 , a trend which is a direct consequence of the 
Maxwellian distribution. When a k = lOOkms" 1 (200kms _1 ), 
the survival probability drops below 10% when gmt ^2.5 AU 
(0.6 AU) for a reference binary mass of 40 M®. 

More input is needed to estimate the net fraction of binaries 
that remain viable BHRP progenitors after the supernova. The 
semimaj or axis and eccen tricity after the supernova are given 
by (e.g.,EiMLMpdsiadlowski 1995) 



/ flSN V 
\flMT / 



1 



A N s+2£w+w 2 



1-A 



NS 



and 



1- 



2 

e SN : 



/ flSN 
\flMT 



(l+£w) 2 + C 2 w 2 



1-A 



(6) 



(7) 



NS 



where w = v^/v/,, and and v^^ are the projections of the 
kick velocity onto, respectively, the directions of the pre- 
supernova orbital velocity and orbital angular momentum. 
After the supernova, the orbit circularizes to a separation of 
flSN.di- = asn(1 "^sn) — assuming conservation of orbital an- 
gular momentum — before the secondary fills its Roche lobe. 
At this stage, we desire the distribution of asN.cir for surviv- 
ing binaries (esN < !)■ To this end, we have adopted the 
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Maxwellian kick distribution, assumed for illustrative pur- 
poses that the amt distribution is logarithmically flat from 
0.2 AU to 10AU, and fixed the total binary mass at 4OM . 
For reasons discussed in § 13.21 this choice for the omt dis- 
tribution may significantly overestimate the number of wide 
binaries. The elementary constraint Osn(1 _e SN) < «mt im- 
plies that flsN.cir < 2a M T for esN < 1; thus, asN.cir < 20 AU in 
this example. Figure |5J the result of a simple Monte Carlo 
integration, shows the cumulative distribution of asN.cir f° r 
systems that remain bound. For = 50kms _1 , lOOkms" 1 , 
and 200 kms" 1 , the fraction of binaries with flsN,cir > 10 AU is 
~2% ~0.4%, and ~0.07%, respecitively. Most systems with 
flSN.dr < 10AU (~2000/?q) will not survive the common- 
envelope phase, as discussed in § 12.21 The cumulative prob- 
ability scales approximately as cr^ 2 5 for asN.cir > 5 AU. An 
increase in the total binary mass from 40M Q to 60M Q in- 
creases the probabilities by at most a factor of 2. In § 13.51 
we combine the results in Fig. [3] with the severe constraints 
of the spiral-in phase in order to more precisely estimate the 
fraction of binaries that survive both the first supernova and 
the spiral-in. 

Supernova mass loss and the neutron-star kick impart some 
velocity to the binary center of mass. Considered indepen- 
dently, the contributions to the systemic speed from mass loss 
and the kick are, respectively, v&Ans(1 + M^/Mf 11 -)' 1 and 
v*(l +M s MT /M NS r 1 (e.g., Appendix B of PRP02). Most bina- 
ries that survive to ultimately form BHRPs will have vt, and 
both less than 100 kms" 1 , since only wide systems can survive 
the common-envelope phase. Therefore, given that Ans < 0-2 
and M^/Mf- T < 0.06, the supernova of the primary's core 
imparts systemic speeds that are typically <10kms _1 . 

3.4. Secondary Winds 

When the secondary's radius is >1000Rq, it may lose most 
of its hydrogen envelope in a wind. A spherically symmetric 
wind that carries the specific orbital angular momentum of 
the star causes the orbit to expand at a rate of d/a = M w /Mb, 
where M w > is the mass loss rate, and Mb is the instanta- 
neous binary mass. During this phase of wind mass loss, the 
Roche-lobe radius of the secondary expands approximately 
in proportion to 85 . If the stellar radius expands with a 
shallower mass dependence, which is quite possible for the 
masses of interest, then it will not be allowed to fill its Roche 
lobe at this stage (PRH03; see also Portegies Zwart, Verbunt, 
& Ergma 1997 and Kalogera & Webbink 1998 for related 
discussions). However, stellar winds from massive stars and 
the effects on binary evolution are very uncertain. The pre- 
scription for the angular momentum taken by the wind can be 
modified so as to limit or reverse the orbital expansion (see 
PRH03), making it possible for the secondary to fill its Roche 
lobe after its mass has been reduced significantly by the wind. 
This may lead to less extreme shrinkage during the common- 
envelope phase (see below), although it seems that the con- 
ditions that permit both Roche-lobe overflow and significant 
mass loss may require fine tuning. This aspect of black-hole 
binary formation deserves considerable attention. 

3.5. Common-Envelope Phase 

As the neutron star spirals in to the stellar envelope, fric- 
tional luminosity is generated at the expense of the orbital 
binding energy. If sufficient energy is available to disperse 
the envelope, the neutron star emerges in a tight orbit about 
the stellar core. A standard energetics argument dWebbink 
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FIG. 4. — Fraction of systems that survive the both the first supernova 
and the spiral-in phase (see text for details), assuming at = 200 km s" 1 . The 
two curves correspond to secondary masses of Mf 1 = 30M Q (solid) and 
50Mq (dashed). The thick dotted line shows a power-law approximation 
intermediate between the two curves curves for 7 > 0.15. 



1984) gives the orbital separation following the spiral-in: 

flCE r L (g SN ) M NS M CtS 
a S N,dr ~ 7 2 MfM„ ' { ) 

where ^sn = A^ns /M™ j <C 1 , M c s and M e >s are the secondary's 
core and envelope mass, and 7 is a combined parameter that 
quantifies the binding energy of the envelope and the effi- 
ciency with which orbital energy is used to expel the enve- 
lope. Over the full range of M^ T , a good approximation to 
eq. (jSJl is 

aw 

— — ~ 0.007 7 (9) 

CSN,cir 

if stellar winds are neglected. Stellar winds prior to Roche- 
lobe overflow reduce the envelope mass but not the core mass. 
A 50% reduction in the envelope mass increases acE/flSN,™- 
by a factor of ~3, allowing a larger number of systems to 
avoid merger. However, as discussed in the last section, a star 
that experiences this degree of wind mass loss is not likely to 
fill its Roche lobe at all. 

The structure parameter 7 varies as a given star evolves. 
Moreover, as the stellar mass is increased, the value of 7 at a 
particular evolutionary stage decreases systematically. Stars 
of mass >25M Q are expected to have 7 < 0.1 in the late 
stages of evolution (Dewi & Tauris 2001; PRH03), and as low 
as ^0.01, giving contraction factors of >1400 from eq. {8J. 
The details of common-envelope evolution are poorly under- 
stood, especially the hydrodynamical aspects. Thus, there is 
considerable uncertainty in the appropriate 7, although it is 
difficult to see how the orbital contraction factor for the sys- 
tems of interest could be much les s tha n 1000 (7 ~ 0.14). 

We now apply the results of § 13.31 and illustrate how the 
fraction of systems that survive the spiral-in depends on 7. 
The core radius is estimated to be (J. Dewi & O. Pols, private 
communication) 

Rc = O.9(M e .. s /lOM Q ) a6 /? . (10) 

We assume a merger results if the core overflows its Roche 
lobe at the end of the spiral-in, i.e., if «ce < ^c,s/'*l(9ce)> 
where qcE = M^s/M c , s . The merger condition then defines 
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the allowed range in asN,cir: 



flSNxir 




We can now extract from Fig. [3] the fraction, fs, of sys- 
tems that survive both the first supernova explosion and the 
common-envelope phase. Results are shown in Fig. 0] for 
Gk = 200 kms" 1 and secondary masses of Mf T = 30Mq and 
50M Q . Smaller values of fs for the 50M Q case are due to the 
larger radius of the exposed core. We have not taken into ac- 
count the proper mass dependence of the supernova survival 
probability. However, this only partially makes up the differ- 
ence between the two curves in Fig.|4] Also shown in Fig.|4] 
is a steep power-law approximation (fs oc 7 2 8 ) to the survival 
fraction that works well for 7 > 0.15, but is clearly an over- 
estimate when 7 < 0.1. If this 'fit' is combined with the 
scaling from § 13.31 we arrive at the simple approximate for- 
mula, 

The survival fraction that results from the above procedure 
assumes that the star can always grow large enough to fill its 
Roche lobe. In fact, the stars of interest probably grow no 
larger than ~2500/?q, which limits the maximum orbital sep- 
aration to ~4OOO/?0 or ~19AU. A smaller maximum radius 
only serves to decrease fs, exacerbating the difficulty in form- 
ing BHRPs. If the orbit contracts by a factor of >1000 dur- 
ing the spiral-in, the resulting orbital separations are <47?©, 
giving orbital periods of <7 hr and relative orbital speeds of 
>700kms _1 for a total binary mass of 10M Q . 

3.6. Supernova of the Secondary and Final Products 

The severity of the orbital contraction in the common- 
envelope phase implies that the orbit must initially be very 
wide and the secondary very evolved in order for the system 
to avoid merger. In most, if not all, cases the exposed core of 
the secondary will have already exhausted its central helium 
supply, so that its remaining lifetime following the spiral-in 
is <10 4 yr. At a mass -loss rate of 10~ 5 -10~ 4 Mq yr -1 (e.g., 
Chiosi & Maed erll986l) . stellar winds only have time to take 
away <1M Q , leading to a <10% widening of the orbit. 

If the eventual core collapse is accompa nied by a suc- 
cessful supernova explosion (see, however, § 12.41 . as much 
as half the mass of the secondary's core may be ejected, 
leading to a large orbital eccentricity. If there is no kick 
to the black hole, the post-supernova eccentricity is esN = 
Abh/(1-Abh), where Abh is the fractional mass loss from 
the binary (eqs. |6) and Q). Black holes probably do re- 
ceive kicks at birth, due to asymmetries in the collapse, explo- 
sion, or in the accretion of fallback material. However, even 
a blac k-hole kick speed of ^100 kms" 1 — as inferr ed for Nova 
Sco (Brandt, Podsiadlowski, & Siaurdsson 1995) — does not 
have a dramatic effect on the orbit, since the pre-supernova 
orbital speeds are >700kms _1 . It is reasonable to suppose 
that, on average, black holes receive a kick momentum that 
is similar to neutron stars, leading to a typical kick speed 
that is perhaps (100-300) x (M N s/M B H)kms _1 , with values 
of ~10-80kms _I for M Bll = 5-10M Q . If typical mass-loss 
and kick parameters are Abh ^ 0.3 and w < 0.2, most newly 
formed BHRPs will have esN S 0-5. Higher eccentricities are 
correlated with larger post-supernova semimajor axes, per- 
haps ~50% larger than the maximum separation following 



the spiral-in phase. In general, we expect that few systems are 
unbound at the time of black-hole formation. 

After the formation of the black hole, further orbital evolu- 
tion is driven by the emission of gravitational radiation. The 
elapsed time before gravitational radiation causes the black 
hole and neutron star to merge is 

r m ~ 2.1 MyrP^M"^ V(l -e 2 ) 7 / 2 , (13) 

where f/, and e are the orbital period (in hours) and eccentric- 
ity, M10 = (Mns+A^bh)/10Mq, and [i\ is the reduced mass in 
Mq. Equation (II 31 is our own fit (good to better than 20% 
for e < 0.9) to the full integral expression for the merger time 
given in iPetersI d!964l) . BHRPs with relatively small eccen- 
tricities and periods of 5-10 hr merge in 10 8 -10 9 yr. 

In general, the mass loss and kick that coincide with the 
formation of the black hole make the largest contribution to 
the systemic velocity of a newly formed BHRP. Because the 
black hole dominates the mass of the binary, most of the kick 
velocity is transferred to the system; an average speed might 
be 20-30kms _I if the kick momenta are similar to neutron 
stars. If the average orbital separation following the spiral-in 
of the neutron star is 2-3^0, then the average pre-supernova 
orbital speed is 800-lOOOkms" 1 . Thus, the mean systemic 
speed due to impulsive mass loss in the second supernova may 
be ~20(ABH/0.1)kms _1 , where Abh is a typical value in this 
context. A typical systemic speed of a few x lOkms" 1 results 
in a vertical scaleheight in the Galactic disk of a few x lOOpc. 
In contrast, it has been suggested that DNSs are born with 
mean systemic speeds of >100kms _1 , acquired mainly after 
the second supernova, giving this population a scaleheight of 
> 1 kpc. 

4. BIRTHRATE ESTIMATE 

We take the total formation rate of massive binaries 
to be the core- collapse supernovae rate, 7?-sn (~10~ 2 yr _1 ; 
CanDe iiaro et al.lll999l) . multiplied by the fraction, fi, ~ 0.5, 
of stars in binaries. Thus, fbHsn closely approximat es th e for- 
mation rate of the primordial binaries discussed in § 13.11 Col- 
lected and summarized below are the efficiencies with which 
the initial ensemble of binaries passes through the various 
stages of evolution described above: 

/mt: The fraction of massive primordial binaries that undergo 
stable mass tra nsfer afte r the primary has left the main se- 
quence (see 13. II and 13.21 for details). For a logarithmi- 
cally flat distribution in initial orbital separations, ~10% of 
the massive binaries satisfy our criteria for stable mass trans- 
fer. Different assumptions that do not depart drastically from 
ours may lead to variations in /mt by a factor of 2 or 3. 

fsH- The fraction of post-transfer secondaries that ultimately 
leave black-hole remnants. In addition to assumptions about 
the distributions of masses and orbital separations for the pri- 
mordial binaries, this fraction depends on the mass threshold 
for black-hole formation and on the amount of matter the sec- 
ondary is able to accept from the primary. Figure suggests 
that /bh ^ 0.04-0.2. 

/s : The fraction of systems that survive the supernova of the 
primary and avoid a merger following the spiral-in phase. Fig- 
ure|4]shows the strong sensitivity of fs to 7 in the specific case 
of o~k = 200 kms" 1 . An approximate scaling of fs with and 
7 is given in eq. (I12> . which overestimates fs for 7 < 0.1. 

Combining these factors, we estimate the Galactic birthrate of 
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Eccentricity 

FIG. 5. — Distribution of orbital periods and eccentricities of newly formed 
BHRPs. The colors were chosen according to the square root of the number 
of systems that enter a given element of a 500 X 500 array. In going from 
magenta to yellow, the number drops by a factor of ~5. 




Eccentricity 

FIG. 6. — Distribution of orbital periods and eccentricities of BHRPs at the 
current epoch. The orbit of each binary has been evolved under the action of 
gravitational radiation losses. The colors are as in Fig. [5] 



BHRPs: 



6 = 5 x 10" 8 yr~ 



K 



fb /ivIT /bH /s 



0.01 yr 1 0.5 0.1 0.1 0.001 



(14) 



Both /mt and /bh could be made larger by a factor of 2 or 
3 over their fiducial values in eq. (114-i . It is then conceivable 
that, if both /mt and /bh are increased significantly, B could 
be as large as ~10~ 6 yr _1 even if /s ~ 10 . However, it is im- 
portant to note that f$ ~ 10~ 3 , which corresponds to 7 > 0.1, 
is quite optimistic in light of recent envelope binding energy 
calculations. We claim that a plausible upper limit to the rate 
is ~10 - ' 7 yr 1 . 

5. GALACTIC POPULATION 

As the final step in our evolutionary study, we have com- 
bined the most important elements of § 3 into a code designed 
to generate a population of newly formed BHRPs. Several 
free parameters are chosen via Monte Carlo methods from 
plausible distributions in such a way as to illustrate the possi- 
ble range of outcomes. Fixed quantities in our simulation in- 
clude the maximum primary mass of 30M Q and the common- 
envelope parameter 7 = 0.1. After a BHRP is formed, we fol- 
low its further orbital evolution due to the emission of gravi- 
tational radiation. Our procedure is outlined below: 

(a) The initial primary mass is chosen from p(M™) oc 
(M™)" 2 - 5 over the range of 8M to 30M Q . The upper limi t is 
intermediate between scenarios I and II discussed in § 13.11 

(b) The secondary mass, Mf B , is drawn from a flat distribu- 
tion in mass ratios. Systems with M™ /M™ < 0.5 are dropped 
since the mass transf er fro m the primary is likely to be dynam- 
ically unstable (see § I3.2> . 

(c) The value of j3 is chosen from the range 0.5-1 with uni- 
form probability. We then compute the secondary mass, Mf T , 
after the phase of stable mass transfer from the primary. We 
adopt 25 Mq as the mass threshold for black-hole formation 
from the secondary, and cut systems with Mf T < 25 Mq. 

(d) Neglecting possible correlations between Mf T and the or- 
bital parameters following the supernova, we assume <7j = 
200km s" 1 and draw flsNck from the cumulative distribution 



shown in Fig. [5] As in Fig. [3] we assume that the distribution 
vanishes at flsN.dr = 20 AU. 

(e) For a given secondary mass, we compute the mass and ra- 
dius of the core (see eqs.Q and 1101 ) that is exposed after the 
common-envelope phase. If the the core overflows its Roche 
lobe following the spiral-in, the system is cut. 

(f) If a system survives the spiral-in, we assume that the sec- 
ondary's core loses 0-2M Q in a stellar wind, chosen with uni- 
form probability. The orbit is widened accordingly. 

(g) We consider both impulsive mass loss and a kick to the 
black hole at the time of core collapse. The fractional mass 
loss is chosen uniformly from the range < Abh < 0.2, and 
the kick is drawn from a Maxwellian distribution with ct^bh = 
200(MNs/MBH)kms _1 , so as to apply an average momentum 
kick that is similar to neutron stars. 

(h) We assume that BHRPs are formed at a constant rate and 
randomly choose a birth time, %, from the range 0-10Gyr. 

(i) After the formation of the BHRP, further orbital evolution 
is driven by gravitation al radi ation, which we calculate using 
the equations of Peters ( 1964). 

We followed 10 9 binaries through the steps above. For each 
successful BHRP that is produced, the event is stored in a 
500 x 500 array representing the distribution of initial orbital 
periods, Pb, and eccentricities, e. Figure[5]shows this distribu- 
tion as a color image, where the colors were chosen according 
to the square root of the number of binaries formed in each el- 
ement of the Pb-e array. In going from magenta to yellow, 
the number per array element decreases by a factor of ~5. 
Most of the newly formed BHRPs have 2hr < Pb < 7hr and 
e < 0.3. The effects of supernova mass loss and black-hole 
kicks lead to a low -probability tail that extends to Pi, ~ 25 hr 
and e ~ 0.8. Not displayed is the distribution of black-hole 
masses, which ranges over Mbh = 5-10M Q , with a mean of 
~7Mq. Our Monte Carlo approach allows us to directly esti- 
mate the birthrate, since we are, in effect, explicitly calculat- 
ing the survival fractions listed in §|4] 

Given the initial Pb, e, and Mbh, the orb i t is ev olved for- 
ward in time using the equations of Peters (1964J); t e is the 
evolutionary time since the birth of the BHRP. The evolution 
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TABLE 1 

Results for Variations in the Standard Model 



Model 


Description 


Birthrate (10 _7 yr') 


Number in Galaxy 


(r m (Myr)> 


<«) 


{^(hr)) 


1 


Standard (7 = 0.1) 


0.63 


5.4 


85 


0.078 


3.5 


2 a ... 


7 = 0.07 


0.03 


0.2 


57 


0.080 


3.2 


3 a ... 


7 = 0.15 


8.80 


103.0 


118 


0.078 


3.8 


4 a ... 


7 = 0.30 


110.00 


3800.0 


345 


0.079 


5.2 


5 b ... 


p(M™) oc CM™)" 2 


0.54 


4.4 


84 


0.078 


3.6 


6 C ... 


p(qm) <x q£ 


0.92 


7.6 


101 


0.079 


3.7 


7 d 


( M MT ),i, - 20Mr^ 


1.50 


26.0 


178 


0.100 


4.1 


8 d ... 


(MMT )th = 3OM 


0.28 


1.3 


47 


0.069 


3.3 


9 e ... 


u k = 100 kms 1 


4.20 


35.0 


84 


0.078 


3.6 


10 c .. 


CT t = 300kms~' 


0.15 


1.3 


85 


0.079 


3.6 


ll f .. 


1.25Rc 


0.14 


2.1 


157 


0.084 


4.9 


12 f .. 


1.5R C 


0.02 


0.4 


280 


0.093 


5.8 


13 s .. 


(AM t .„,)wmd=0 


0.63 


2.5 


40 


0.072 


3.0 


14 h .. 


A BH = 


0.63 


2.8 


44 


0.042 


2.9 


15'... 


f,t,BH = km s" 1 


0.63 


4.8 


76 


0.058 


3.6 


16'... 


o-t.BH = 400(Mns / Mbh ) km s -1 


0.63 


5.9 


93 


0.116 


3.7 



"Envelope binding-energy parameter associated with the spiral-in phase (eq.l8*l). 

b The initial mass function for the primary star. 

c The distribution of mass ratios for the primordial binaries. 

d Minimum post-transfer secondary mass that leaves a black-hole remnant. 

e The one-dimensional velocity dispersion for neutron-star kicks. 

The radius of the secondary's core, nominally given by ea.fiol. was arbitrarily increased by factors of 1.25 and 1.5 to allow for the effects of 
thermal bloating. 

g Maximum mass loss via a stellar wind from the exposed core of the secondary. 
h Fractional mass lost from the binary during black hole formation. 
'The one-dimensional velocity dispersion for black-hole kicks.. 



is terminated when the system merges for tb+t e < lOGyr, or 
when tt,+t e = lOGyr is reached. For each BHRP that traverses 
the Pb-e plane, we store the amount of time spent in each el- 
ement of the 500 x 500 array; these times are accumulated 
for all the systems. The result is an array of total dwell times 
that represent the probability density of finding a system at 
the current epoch with a given Pb and e. Figure[6]displays the 
expected current-epoch distribution of binary parameters. 

Summing the dwell times over the entire array, dividing 
by the total number of BHRPs formed successfully from 10 9 
primordial binaries, and multiplying by the BHRP birthrate, 
we estimate the number of BHRPs that currently inhabit the 
Galaxy. In our standard model listed in the steps above, we 
find 6~6x 10" 8 yr _1 and estimate that ~5 BHRPs populate 
the Galaxy. This estimate can also be obtained upon exam- 
ination of Fig. [6] where it is evident that the average orbital 
period is 3-4 hr, corresponding to a short merger time of ~(4— 
8) x 10 7 yr. 

Table 1 gives an indication of how our results vary when 
we consider a range of plausible input assumptions that differ 
from those in our standard model. For each quantity or func- 
tion that is varied, all else is left as in the standard model. As 
already demonstrated, the BHRP birthrate and total number 
in the Galaxy depend most strongly on the common-envelope 
parameter 7. In the very extreme case of 7 = 0.3, the birthrate 
and total number are ^10~ 5 yr" 1 and ~4000, respectively. For 
7 = 0.15, which is on the upper edge of plausibility, the total 
number is ~100. The column labeled (r m (Myr)}, an average 
merger time, is defined to give the total number when multi- 
plied by the birthrate in column 3. In columns 6 and 7, (e) 
and (P(hr)) are traditional means of the distribution. 

6. PULSAR DETECTION AND TIMING 



Our results suggest that at most ~100, but probably <10, 
BHRPs currently populate the Galactic disk. We now con- 
sider what fraction of these systems are detectable in radio 
pulsar surveys. A plausible first guess is that, as a class, 
BHRPs are similar observationally to their DNS cousins. A 
number of papers (e.g., Kalogera et al. 2001 and references 
therein) indicate that < 1 % of the total Galactic DNS popula- 
tion has been discovered among several extensive pulsar sur- 
veys conducted over the past 30 yr. This fraction incorporates 
the sensitivity and sky coverage of each survey, the expected 
spatial and luminosity distributions of DNSs, and the fraction 
(~20-30%) of pulsars beamed toward the Earth. Such a small 
detection efficiency is consistent with not yet having discov- 
ered a BHRP. If indeed there are ~100 BHRPs in the Galaxy, 
then ~1 system may be within reach of ongoing surveys us- 
ing the Parkes and Green Bank radio telescopes. However, 
we view this as being very optimistic, since 1% is probably 
a firm upper limit to the detection efficiency with current ra- 
dio telescopes, and because we predict a substantially smaller 
number of BHRPs. Realistically, there are likely to be impor- 
tant differences between the BHRP and DNS populations that 
impact the BHRP detection likelihood, three of which we now 
discuss in turn. 

The recycling histories of B HRPs and DNSs may be dis- 
tinct. As discussed in § 12.21 it is possible that pulsars in 
BHRPs have longer spin periods, P s , and perhaps also larger 
surface magnetic fields, B. The spin-down luminosity scales 
as B 2 /P*. Systematically longer spin periods may be com- 
pensated for by larger magnetic fields, giving BHRPs lumi- 
nosities similar to DNSs. 

For reasons given in § 13.61 the mean systemic speed (e.g., 
~30kms _1 ) of newly formed BHRPs may be smaller than 
for DNSs (>100kms _1 ), resulting in a smaller scaleheight of 
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BHRPs in the Galactic disk: <500pc for BHRPs, compared 
to >1 kpc forDNSs. Overall, we expect that the smaller scale- 
height of BHRPs somewhat increases their discovery proba- 
bility relative to DNSs, but not by a large factor. It is also 
important to note that if BHRPs are more concentrated to- 
ward the Galactic midplane, pulse smearing due to dispersion 
in the interstellar medium, which decreases the sensitivity, is 
more significant than for DNSs. However, the modest relative 
reduction in sensitivity is partly canceled if pulsars in BHRPs 
have longer spin periods. 

Perhaps the most robust difference between BHRPs and 
DNSs, in regard to their detectability, is that the pulsars in 
BHRPs should, on average, have larger orbital accelerations, 
because of their more massive black-hole companions and 
generally short orbital periods. Acceleration of the pulsar 
causes the measured pulse frequency to drift over the integra- 
tion time of the observation, thus smearing the pulsar's signal 
in the Fourier p ower spectrum of the data and reducing it s 
amplitude (e.g., Anderson 1993; Johnston & Kulkarni 1991). 
For black-hole masses in the range of 5-10M Q , the acceler- 
ation of a pulsar in a BHRP is ~2 times that in a DNS with 
the same orbital period and binary inclination, assuming that 
the orbits have small eccentricities. A search that explicitly 
probes for accelerated pulsars may be required in order to de- 
tect a BHRP. Such acceleration searches are computationally 
expensive, but are now done routinely. 

Once a BHRP is discovered, the rewarding task of timing 
the pulsar begins. The initial Keplerian orbital parameters will 
be readily determined with high precision. Since BHRP orbits 
are probably very compact, the effects of relativistic gravity 
will be quite pronounced. Relativistic effects that have been 
measured in binary pulsars include orbital precession, the 
combined effect of the second-order Doppler shift and grav- 
itational redshift (Einstein delay), the Shapiro propagation 
delay, orbital contraction due to gravity-wave emission, and 
the infl uence of geodetic precess i on on the measured pulse 
profile JDamour & Tavlorl 1 19921 FWeisber glTTavloil 120021: 
IStairs et alJl2004l) . Each of these effects is a priori more sig- 
nificant in a BHRP than in a DNS, as are d elays due to rel- 
ativis tic aberration and light bending (e.g., Wex & Kopeikin 
1999). The Keplerian parameters, coupled with the measure- 
ment of the Einstein or Shaprio delay, will yield the mass of 
the black hole, which is fundamentally important for study- 
ing models of black-hole formation, in particular because the 



black holes in BHRPs have not accreted mass. A BHRP also 
offers the tantalizing possibility of measuring strong-field ef- 
fects of relativistic gravity and providing stringent constraints 
on theories alternative to Einstein's general relavity. The high 
sensitivity of the planned Square Kilometer Array may be re- 
quired to tap into these more exotic problems (e.g., Kramer et 
al. 2004 and references therein) — and perhaps to provide the 
first BHRP detection. 

7. CONCLUSIONS 

We would like to underscore three main conclusions of our 
study. (1) There are major uncertainties in the formation of 
BHRPs, linked to important details in the life and death of 
stars more massive than ~20M Q , but mainly to a poor un- 
derstanding of the common-envelope phase. Under quite op- 
timistic circumstances, we expect that the Galactic BHRP 
birthrate is ^10~ 7 yr _1 . The birthrate may be much lower, 
and possibly zero. The effects of dramatic orbital shrinkage 
during the common-envelope phase, prodigious stellar winds, 
and hypercritical acc retion may preclude BHRP formation 
altogether (see § 12. 2\ . (2) Plausible, but possibly minimal, 
common-envelope shrinkage factors of ~1000 lead to initial 
BHRP orbital periods of 2-10 hr and an average merger time 
of ^10 8 yr. (3) The short merger times and low birthrate im- 
ply that there are most likely no more than ~10 BHRPs cur- 
rently in the Galactic disk, making the prospects for their de- 
tection somewhat bleak. As a direct comparison, we estimate 
that there is probably no more than 1 BHRP in the Galaxy 
for every 100-1000 double neutron stars, of which only 8 are 
currently known. 

Our simulations suggest that most BHRPs at the current 
epoch have orbital periods of 1-6 hr. Due to supernova mass 
loss and possible black-hole kicks, a typical eccentricity may 
be ~0. 1 . The acceleration of the pulsar in such a binary is 
large, but does not greatly hinder its detection given current 
computational resources. If a BHRP is discovered, timing 
measurements would provide a wealth of new constraints on 
theories of gravitation. 

We would like to thank J. Dewi, R. Di Stefano, A. Loeb, 
and S. Ransom for valuable discussions. EP was supported 
by NASA and the Chandra Postdoctoral Fellowship program 
through grant number PF2-30024. 



REFERENCES 



Alpar, M. A., Cheng, A. R, Ruderman, M. A., & Shaham, J. 1982, Nature, 
300, 728 

Anderson, S. B. 1993, Ph.D. thesis, Caltech 

Arzoumanian, Z., Cordes, J. M., & Wasserman, I. 1999, ApJ, 520, 696 
Arzoumanian, Z., Chernoff, D. F., & Cordes, J. M. 2002, ApJ, 568, 289 
Belczynski, K., Kalogera, V., & Bulik, T. 2002, ApJ, 572, 407 
Bhattacharya, D. & van den Heuvel, E. P. J. 1991, Phys. Rep., 203, 1 
Brandt, N. & Podsiadlowski, P. 1995, MNRAS, 274, 461 
Brandt, W. N., Podsiadlowski, P., & Sigurdsson, S. 1995, MNRAS, 277, L35 
Brown, G. E., 1995, ApJ, 440, 270 

Brown, G. E., Heger, A., Langer, N, Lee, C.-H., Wellstein, S., & Bethe, H. 

2001, New Astronomy, 6, 457 
Burgay, M., et al. 2003, Nature, 426, 531 
Cappellaro, E., Evans, R., & Turatto, M. 1999, A&A, 351, 459 
Chevalier, R. A. 1993, ApJ, 411, L33 
Chevalier, R. A. 1996, ApJ, 459, 322 
Chiosi, C, & Maeder, A. 1986, ARA&A, 24, 329 
Damour, T., & Taylor, J. H. 1992, Phys. Rev. D, 45, 1840 
Dewi, J. D. M. & Tauris, T. M. 2001, ASPConf. Ser. 229, Evolution of Binary 

and Multiple Star Systems, ed. Ph. Podsiadlowski, S. Rappaport, A. R. 

King, F. D'Antona, & L. Burderi (San Francisco: ASP), 255 



Dewi, J. D. M. & Pols, O. R. 2003, MNRAS, 344, 629 

Dymnikova, I. G. 1986, IAU Symp. 114, Relativity in Celestial Mechanics 
and Astrometry: High Precision Dynamical Theories and Observational 
Verifications, ed. J. Kovalevsky & V. A. Brumberg (Dordrecht: Reidel), 
411 

Eggleton, P. P. 1983, ApJ, 268, 368 
Faulkner, A. J., et al. 2005, ApJ, 618, LI 19 
Fryer, C. L. & Kalogera, V. 2001, ApJ, 554, 548 

Goicoechea, L. J., Mediavilla, E., Buitrago, J., & Atrio, F. 1992, MNRAS, 
259, 281 

Gourgoulhon, E., & Haensel, P. 1993, A&A, 271, 187 
Hansen, B. M. S. & Phinney, E. S. 1997, MNRAS, 291, 569 
Hulse, R. A. & Taylor, J. H. 1975, ApJ, 195, L51 

Humphreys, R. M. 1 984, in Observational Test of Stellar Evolution Theory, 

ed. A. Maeder & A. Renzini (Dordrecht, Reidel), 279 
Hurley, J. R„ Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 
Ivanova, N., Belczynski, K, Kalogera, V., Rasio, F. A., & Taam, R. E. 2003, 

ApJ, 592, 475 
Johnston, H. M. & Kulkarni, S. R. 1991, ApJ, 368, 504 
Jonker, P. G., & Nelemans, G. 2004, MNRAS, 354, 355 
Joss, P. C. & Rappaport, S. A. 1983, Nature, 304, 419 



10 



Pfahl, Podsiadlowski, & Rappaport 



Kalogera, V. & Webbink, R. F. 1998, ApJ, 493, 351 

Kalogera, V., Narayan, R., Spergel, D. N., & Taylor, J. H. 2001, ApJ, 556, 
340 

Kalogera, V., et al. 2004, ApJ, 601, LI 79 

Kramer, M., Backer, D. C, Cordes, J. M., Lazio, T. J. W., Stappers, B. W., & 

Johnston, S. 2004, New Astronomy Review, 48, 993 
Laguna, P. & Wolszczan, A. 1997, ApJ, 486, L27 
Langer N., & Maeder A., 1995, A&A, 295, 685 

Lipunov, V. M., Postnov, K. A., Prokhorov, M. E., & Osminkin, E. Y. 1994, 

ApJ, 423, L121 
Lyne, A. G., et al. 2004, Science, 303, 1 153 
Narayan, R., Piran, T., & Shemi, A. 1991, ApJ, 379, L17 
Peters, P. C. 1964, Physical Review, 136, 1224 

Pfahl, E., Rappaport, S., & Podsiadlowski, P. 2002, ApJ, 573, 283 (PRP02) 
Pfahl, E., Rappaport, S., Podsiadlowski, P., & Spruit, H. 2002, ApJ, 574, 364 
Podsiadlowski, Ph., Nomoto, K., Maeda, K., Nakamura, T., Mazzali, P., & 

Schmidt, B. 2002, ApJ, 567, 491 
Podsiadlowski, P., Rappaport, S., & Han, Z. 2003, MNRAS, 341, 385 

(PRH03) 

Podsiadlowski, Ph., Langer, N, Poelarends, A., Rappaport, S., Heger, A., & 
Pfahl, E. 2004, ApJ, 612, 104 



Portegies Zwart, S. F., Verbunt, E, & Ergma, E. 1997, A&A, 321, 207 
Ransom, S. M., et al. 2004, ApJ, 609, L71 

Sigurdsson, S. 2003, in ASP Conf. Ser. 302, Radio Pulsars, ed. M. Bailes, D. 

J. Nice, & S. E. Thorsett (San Francisco: ASP), 391 
Sipior, M. S., & Sigurdsson, S. 2002, ApJ, 572, 962 

Sipior, M. S., Portegies Zwart, S., & Nelemans, G. 2004, MNRAS, submitted, 

astro-ph/0407268 
Smarr, L. L. & Blandford, R. 1976, ApJ, 207, 574 
Stairs, I. H. 2004, Science, 304, 547 

Stairs, I. H., Thorsett, S. E., & Arzoumanian, Z. 2004, Phys. Rev. Lett., 93, 
141101 

Tauris, T. M., & van den Heuvel, E. P. J. 2003, to appear in "Compact Stellar 
X-ray Sources, eds. W. H. G. Lewin & M. van der Klis, astro-ph/0303456 
Voss, R. & Tauris, T. M. 2003, MNRAS, 342, 1 169 
Webbink, R. F. 1984, ApJ, 277, 355 
Weisberg, J. M., & Taylor, J. H. 2002, ApJ, 576, 942 
Wellstein, S. & Langer, N. 1999, A&A, 350, 148 
Wex, N. & Kopeikin, S. M. 1999, ApJ, 514, 388 
Woosely, S. E., Weaver, T. A. 1995, ApJS, 101, 181 



